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ABSTRACT 

Many small bodies in the solar system are believed to be rubble piles, a collec- 
tion of smaller elements separated by voids. We propose a model for the structure 
of a self-gravitating rubble pile. Static friction prevents its elements from sliding 
relative to each other. Stresses are concentrated around points of contact between 
individual elements. The effective dimensionless rigidity, fi ru bbie, is related to that 
of a monolithic body of similar composition and size, fi by fi r ubUe ~ /t^ey 1 ^ \ 
where ey ~ 10 -2 is the yield strain. This represents a reduction in effective 
rigidity below the maximum radius, R m ax ~ [fJ-cy / \Gp 2 )] 1 / 2 ~ 10 3 km, at which 
a rubble pile can exist. Densities derived for binary near-Earth asteroids im- 
ply that they are rubble piles. As a consequence, their tidal evolution proceeds 
10 3 to 10 4 times faster than it would if they were monoliths. This accounts for 
both the sizes of their semimajor axes and their small orbital eccentricities. We 
show that our model for the rigidity of rubble piles is compatible with laboratory 
experiment in sand. 



Subject headings: asteroids 



1. Introduction 



Rubble piles are bodies composed of smaller elements separated by voids. There is 
compelling evidence that at least some small solar system bodies are rubble piles bound by 
gravity. Their telltale signature is a mean density below that of their constituent elements. 
Examples include: four icy satellites of Satur n, the coorbitals, Janus and Epimetheus, an d the 
F-ring shepherds, Prometheus and Pandora (IJacobson fc French! 2004l; |Porco et al.ll2007l). the 
rocky main belt asteroid s, C-type 253 Mathilde ( Veverka et al. 1997 ; lYeomans et al. 1997 ) 
and M -type 22 Kalliope (IMargot fc Brown! 120031 ). and the binary near-Earth asteroid 1999 
KW4 (lOstro et al.l 120061 ) . The largest of these bodies, Janus and Kalliope, have dimensions 



-2- 



of order 100 km. It is unclear whether larger rubble piles exist or whether all smaller bodies 
are rubble piles. 

Intuitively we expect a rubble pile to be weaker than a monolithic body of the same 
composition. Thus tidal dissipation at a rate that is more rapid than typical for a monolith 



is considered evidence for a rubble pile (IMargot fc Brown! 120031 ). The orbits of binary near 



Earth asteroids are prime examples; the sizes of their semimajor axes and their low orbital 
eccentricities suggest that they are evolving tidally at rates between 10 3 and 10 4 times faster 



than e stimates for monolithic bodies of similar size would predict ( jWalsh fc Richardson 



2006al ). In what follows, we provide a theoretical basis for estimating tidal dissipation rates 



in rubble piles and show that it can account for this large speedup of tidal evolution. 

The plan of our paper is as follows. In §21 we formulate a quantitative theory for the 
effective rigidity of a self-gravitating rubble pile and demonstrate that it is due to voids 
rather than cracks. Limits on the sizes of rubble piles are derived in |j3l In we apply our 
theory to the tidal evolution of binary near-Earth asteroids. 



2. Effective Elastic Modulus Of A Rubble Pile 

We begin by reviewing the tidal response of a uniform body of density p, rigidity p, and 
radius R. As is customary, we define the dimensionless rigidity by p; 

2gpR 

Next we show that p is the ratio of the fluid strain to the elastic strainQ 

We assume that the tidal force, /, is weak in comparison to the cohesive force of the 
body's self gravity, gM, where g ~ GpR. If the body were fluid, p = 0, it would suffer a 
strain 

£ '~^- (2) 

whereas if it were elastic but lacked self-gravity, g = 0, the strain would be 

To order of magnitude, the ratio between expressions <^ and (|3J) reproduces p given by 
equation (Tj[|). 



Arguments in this section are order of magnitude only. 
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How does the tidal response of a rubble pile differ from that of a monolith? To answer 
this question, we investigate some simple models. 



2.1. cracks do not matter 

Normal stresses are seamlessly transmitted across cracks, so a body's response to weak 
tides is unaffected by cracks. Consider a body of radius R composed of cubical elements 
whose sides have length r <C R. The ratio of the weight of a single cube, gM(r/R) 3 , to the 
divergence of tidal stress acting on its volume, f(r/R) 3 , is just gM/ f. Thus a coefficient 
of static friction larger than f/gM would suffice to prevent the cubes from slipping relative 
to each other. Coefficients of static friction for rocks and dry ice are of order unity, and 
f/gM ~ (R/a) 3 for an equal mass binary with semi-major axis a. 



2.2. voids are key 

2.2.1. uniform spheres 

Next we consider a body of radius R composed of identical spheres of radius r <C R. 
Its mean density p ~ 0.7p. A typical cross section cuts (R/r) 2 small spheres each of which 
transmits forces F(r/R) 2 to its neighbors, where F ~ gM + / includes both tidal forces 
and self gravity. In so doing, each small sphere undergoes a linear distortion d~x and forms 
contact surfaces with its neighbors of area Sxr. Within (Sxr) 1 ^ 2 of each contact, the strain 
is of order (Sx/r) 1 ^ 2 so 

F (jA* ~ i*r 1 ' 2 6!i?' 2 . (4) 
The average strain is just Sx/r, where from equation (j3J) 



5x ( F \ 



2/3 



(5) 



Most of this strain is due to the body's self-gravity. To isolate the tidal strain, we expand 
F 2 ! 3 in equation (jSJ) around F ~ gM to obtain 

1/3 

(6) 



f ( P 



fxR 2 \gpR 

Thus the effective dimensionless tidal rigidity of a body composed of identical spheres is 

^spheres ~ \ T^J ~ ' ^ 
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This result is equivalent to that originally established by iDuffy fc Mindlinl (119571 ) 



2. 2. 2. irregular fragments 

Natural rubble piles are likely to be composed of irregularly shaped elements whose 
surfaces have local radii of curvature, f , that are much smaller than the elements' sizes, 
r. Compared to rubble piles composes of spheres, contact areas would be reduced, stress 
concentrations increased, and the effective rigidity lowered. A simple modification of the 
derivation given in 12.2.11 suffices to evaluate the effective rigidity of a rubble pile, \x ru bUe- 
Each element still transmits its share of the total force. However, f must replace r on the 
right hand side of equation (j4]). Thus now 

5x ( F \ 2/3 /r\i/3 ^ 
r \fj,R 2 J \fJ 

Continuing as before, we find that the average strain across the rubble pile is increased by 
the factor (r/f) 1 / 3 with the consequence that the effective rigidity now reads 

/f\ 1/3 2 3 /r\ 1/3 

firubble ~ ^spheres f ~ J ~ ( ~ J (9) 

The sharper the contact points, the softer the rubble pile, up to a limit at which the stress 
near the contact surfaces reaches the material's yield stress ay, or yield strain e Y = oy/jti. 
This limit is met at 



f 1 



(10) 



m) 

Sharper contact points than allowed by equation ( flOl) would be dulled by material flow or 
failure. Therefore, 

/ LL \ V2 

^rubble > — • (11) 



Experimentally it is generally found that the effective rigidity of a granular mate rial scales 



in direct proportion to the square root of the confining pressure. iGoddardl (119901 ) provides 
an explanation for this scaling which is similar to ours. 

Equations (|7j) and (Q demonstrate that the effective rigidity of a rubble pile is smaller 
than that of a monolithic body of the same size. The reduction in rigidity is independent of 
the sizes of the elements into which the body is divided. It arises from the concentration of 
stresses due to the presence of voids. 
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TTT 



/i 



19 fJL 

no sliding 





^rubble ~ (A/ey) 1 ^ 2 
< 1000 km 



Fig. 1. — Three simple models of fragmented bodies. Upper row depicts a body composed 
of cubical elements. There are no voids. Static friction prevents the elements from sliding 
relative to each other. Its effective rigidity is identical to that of a monolith. Middle row 
shows a body composed of spherical elements. Voids are present. Stresses concentrate near 
contacts between elements. The effective rigidity is smaller than that of a monolith and is 
independent of the sizes of the elements. Such an idealized configuration requires the spheres 
to be sufficiently smooth. If made of rock or ice, its radius could not be larger than about 
10 km. Bottom row illustrates a more realistic rubble pile composed of irregular elements. 
Sharper contact points increase stress concentration more than for a body composed of 
spherical elements. Accordingly its effective rigidity is further decreased. Radii of rubble 
piles composed of rock or ice cannot be larger than about 1000 km. 
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2.3. energy considerations 

We re-derive equation ([9]) based on energy considerations. Strains of order (Sx/f) 1 ^ 2 are 
attained within a distance (Sxf) 1 ^ 2 of the contacts between individual elements. Thus the 
elastic energy stored within the rubble pile satisfies 

SE ~ ^{Sxf'H 1 / 2 (jj 3 . (12) 

We can also express SE in terms of the effective dimensionless rigidity, jX TU bbie and the average 
strain in the rubble pile, Sx/r as 

/S x \ 2 

6E~fi rubble gpRl—J R 3 . (13) 
Equating the expressions for SE given in equations (TT2l) and (1131) . we arrive at 



„8x l/2 r 1 / 2 

^rubble ~ H • (14) 

r 



Finally, by using equation (JSj) to eliminate Sx, we recover equation (jUJ) 



3. Critical Sizes For Rubble Piles 
3.1. mechanical limits 

At 

R = R,~(0'\ (!5) 

which corresponds to /2 ~ e^ 3 , fi rv hble ~ ^spheres ~ £y 2 and f/r ~ 1. For nominal values of 
AWfc ~ 5 x 10 11 dyne cm -2 , // ice «3x 10 10 dyne cm -2 , ey ~ 10" 2 , i?* ~ 10 km for rubble 
piles composed of either rock or ice. Moreover, jX ru bbie ~ 10 4 as compared to /t ~ 10 6 for a 
monolith of radius R*. At 

/ \l/2 
-Rmax = ( ^2£ J ' 

which corresponds to jl ~ e y x and f/r ~ e^ 2 , the contact areas are comparable to the surface 
areas of individual elements so fi ru bbie ~ A- With nominal parameters, -R max ~ 10 3 km and 

/^rufcftZe ~ At ~ 10 2 . 

A body with R < R* would avoid elastic failure if it were composed of identical spheres. 
For elastic failure would occur at points of contact among spheres. More generally, 
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we would expect the voids in rubble piles to occupy a decreasing fraction of the volume with 
increasing R up to R = R max - At R ma x, the average interior pressure gpR max ~ ay, so voids 
could only exist near the surface. 



3.2. thermal limits 

Rubble piles should be more common among smaller bodies because they cool more 
rapidly than larger ones and therefore are less likely to have been melted. Thermal diffusiv- 
ities of rock and ice are of order 10~ 2 cm 2 s _1 , which implies 

*-~ 3x1010 (lo^) 2y - (17) 

Even bodies as small as R* ~ 10 km might have been melted if they formed early and 
were endowed with short lived radioactive nuclides. On the other hand, bodies as large as 
Rmax ~ 3 x 10 2 km which were fragmented by collisions after the short lived radioactive 
nuclides had decayed could have avoided melting. 



4. Implications For Tidal Evolution 



Tides play crucial role in orbital and spin evolution of binaries. Here we focus on the 
evolution after the secondary's spin has become synchronous with the mean orbital angular 
velocity while the primary's spin remains much faster. In this case, tides raised on the 
primary cause both the semimajor axis, a, and the orbital eccentricity, e, to grow while 
those raised on the se condary have negligible effect on the sem imajor axis and cause the 
eccentricity to decay (iGoldreichl 1 19631 ; iGoldreich fc Soterl Il966l ) . Relevant expressions for 
e«l are: 

1 da 
a dt 

and 

1 de 
e dt 

for tides raised on the primary, and 



Q P M P 



57 k p M s 
~8Q~ P M~ P 



R, 



R„ 



n 



n . 



(19) 



lde _ 21 k s M p ( R s 
edt " ~T<^Ms V^" 



n . 



(20) 
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for tides raised on the secondarvn 



Tidal evolution rates depend on two parameters, k and Q. The estimation of k for 
monoliths involves little uncertainty. For a body of uniform density, 

1.5 



k 



The estimation of Q is more uncertain. Available evidence suggests that Q 
monolithic bodies (jGoldreich fc Soterlll966l ). 



(21) 
10 2 for 



4.1. Semimajor axis evolution in binary rubble piles 

Semimajor axis evolution is driven by the transfer of angular momentum from the spin 
of the primary to the orbit. Below, we estimate timescales for the semimajor axes of some 
well observed binary NEAs to have evolved from much smaller initial values to their current 
ones. Integrating equation ( 1181) . we obtain 

m 2 Q v M„ ( a \ 5 1 



39 k p M s \R P J n 

We compare timescales for models in which the bodies are assumed to be monolithic solids, 
fluids, and rubble piles. We set Q = 100 in each case. 

As the entries in table [TJ demonstrate, the timescale for semimajor axis evolution is 
measured in Gyrs for monoliths, years for fluids, and Myrs for rubble piles. Only the latter 



is con sistent with estimates of 10 Myr for the dynamical life time of NEAs ( iGladman et al. 



20001 ) . Since it is plausible that the stress concentration in rubble piles results in Q < 100, 



the ages we estimate for rubble piles should be viewed as upper limits. 



4.2. Comparison with experiments in sand 

The effective rigidity of our model rubble pile, fj, ru bbie, is proportional to the square root 
of the confining pressure and independent of the size of the individual elements. Laboratory 



2 Subscripts p and s d enote primary and second ary. We adopt standard notations for tidal Love number, 
k, and quality factor, Q (jMurrav fc Dermottll200(I K Q" 1 is a stand-in for sin2<5, where S is the tidal phase 
lag. 
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measurements of the shear velocity, c s = a/ p/p, in sand as a function of pr essure provide 



a useful calibration. The data on c s (p) plotted in figure 1 of iGoddardl (119901 ) are replotted 
in our figure 2. On the figure's upper boundary we display the radius of an asteroid whose 
average internal pressure 

/ R \ 2 

P = (47r/15)Gp 2 i? 2 S 2.2 x 10 3 I —— J dyne cm" 2 (23) 

with p = 2 g cm -3 corresponds to that given on the lower boundary. The range of pressures 
covered in the experiments on sand correspond to those inside asteroids with radii from 
10-40 km. The right-hand boundary of the figure shows the effective rigidity corresponding 
to the shear velocity. It is well-fit by the solid line which is derived from our expression for 
effective rigidity with ey = 0.17. This should not be taken as evidence that the yield stain 
of sand is 0.17 since our formula is only accurate to order of magnitude. However, it does 
suggest that the ages we estimate in table [Q may be a factor of a few too large. The dashed 
line indicates the higher effective rigidity of a body composed of uniform quartz spheres. 

Next we compare data on the rigidity of sand with that on the effective rigidity of 
NEAs. To do so, we assume that the semimajor axes of binary NEAs have evolved from 
much smaller initial values over ~ 1 Myr with a tidal Q = 100. Then we use equations (pQ), 
( I18p . and (j2ip to evaluate the effective tidal rigidity of the primary for each of the binaries in 
table CD These rigidities are plotted as x's on figure Although the scatter is large, probably 
dominated by our assumption of a uniform age, the data fit nicely on the extrapolation to 
low pressure of the data from the experiments on sand. 



4.3. rates of eccentricity evolution in binary rubble piles 



Binary near-Earth asteroids are thou ght to form by YorrJ 3 



ruption and consequently be rubble piles (jWalsh fc Richardson 



spin u p and/or tidal dis- 
2006bl ). Most have nearly 



circular orbits from which IWalsh fc Richardson! (j2006al ) argue that tidal damping of their 
orbital eccentricities proceeds 3 to 4 orders of magnitude faster than would be expected for 
binary monoliths of comparable size. A significant fraction of this increase must be due to 
the reduced rigidity of a rubble pile as compared to a monolith since Q cannot be smaller 
than unity. 



3 Y orp spin up has been measured for near-Earth asteroid 2000 PH5 (jLowrv et al. 
20071) . 



2007: 



Taylor et al 
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Effective Radius (km) 
0.1 0.3 0.5 1 2 5 10 20 40 




Pressure (dyn/cm 2 ) 



Fig. 2. — Comparison of our mod el for the e ffecti ve rigidity of rubble piles with that from 
experiments on sand taken from iGoddardl (jl990l ). Shear wave velocity as a function of 
pressure in saturated, dry and drained Ottawa sands is shown by open circles, solid circles, 
and triangles, respectively. Effective rigidities of NEAs, inferred by assuming binary ages of 
lMyr and Q = 100, are plotted against the primary diameter and marked by x's. 
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For rocky bodies, scaling from the tidal Love number of the Moon, kuoon ~ 0.03, 

f kin "\ ^ 

p, « 1.5 x 10 8 -5- , (24) 



v R t 

which corresponds to /i ~ 5 x 10 11 dyne cm -2 . Thus from equation (II ip with ey = 10 -2 , we 
obtain 

J- < (^) 1/2 « 10'^. (25) 

prubble 

Since typical secondaries among near earth asteroid binaries have radii of a few tenths of a 
kilometer, much if not all of the increase in the inferred rates of eccentricity damping might 
be due to an increase of k. However, it would not be surprising if a contribution came from 
a reduction of Q. 

We note that close encounters with Earth or other planets might reset the eccentricities 
of binary NEAs on timescales comparable to those at which they evolve under tides. This 
issue deserves investigation. 



4-3.1. conditions for eccentricity damping in binary asteroids 



If both primary and secondary were strength rather than gravity dominated (/j > 1) 
then the ratio of the rates of eccentricity excitation and damping would be 

19 f p s \ 2 R s jx s Q s 



^ = ^i^i^i. (26) 

28 \PpJ R p p p Q p 

For monoliths of identical composition, this ratio reduces to 

'^monolith _ Q r-> /O ' \ ) 

lo K s LJp 

Thus for Qs/Qp = 10 eccentricity damping would require R p /R s < 1.47 corresponding to a 
mass ratio less than 3.20. For primary and secondary composed of spherical elements with 
identical compositions and Q's, the ratio reads 

19 [R P \ l/3 

TZ-spheres — 28 ( J ' 



Identical compositions do not guarantee identical Qs, because the latter may also depend on strain, 
strain rate, temperature, and pressure. 
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so eccentricity would damp for R p /R s < 3.2 corresponding to a mass ratio below 33. Finally, 
for rubble piles composed of irregular elements of identical compositions and Q's, 

19 

Ttrubble = 7777 (29) 
Zo 

so eccentricity would damp for all mass ratios. 

It is clear that eccentricity damping is more likely for binary rubble piles than for binary 
monoliths especially when the mass ratio is not large. However, given the uncertainties in the 
relative values of the primary's and secondary's ji ru bbie an d Q, eccentricity growth remains a 
possibility, in particular for large mass ratios. 

This research was supported in part by an NSF grant and a NASA grant. RS is an Alfred 
P. Sloan Fellow, and a Packard Fellow. We thank Hiroo Kanamori for valuable advice. 
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Asteroid 


Orbital 


Semimajor 


Primary 


Secondary 


Monolith 


Fluid 


Rubble Pile 


name 


period (days) 


axis (km) 


diameter (km) 


diameter (km) 


age (Gyr) 


age (yr) 


age (Myr) 


(66391) 1999 KW 4 


0.73 


2.5 


1.2 


0.4 


15 


37 


7.5 


1999 HFi 


0.58 


7.0 


3.5 


0.8 


3.6 


74 


5.2 


(5381) Sekhmct 


0.52 


1.5 


1.0 


0.3 


4.2 


7.0 


1.7 


(66063) 1998 ROi 


0.6 


1.4 


0.8 


0.38 


4.1 


4.4 


1.3 


1996 FG 3 


0.67 


2.6 


1.5 


0.47 


4.3 


16 


2.7 


(88710) 2001 SL 9 


0.68 


1.4 


0.8 


0.22 


24 


26 


7.8 


1994 AWi 


0.93 


2.3 


1 


0.5 


14 


23 


5.6 


2003 YTi 


1.2 


2.7 


1 


0.18 


880 


1500 


360 


(35107) 1991 VH 


1.4 


3.2 


1.2 


0.44 


74 


180 


36 


2000 DPiot 


1.8 


2.6 


0.8 


0.3 


540 


580 


180 


(65803) Didymos 


0.49 


1.1 


0.8 


0.17 


11 


12 


3.7 


(5407) 1992 AX 


0.56 


6.8 


3.9 


0.78 


2.1 


54 


3.4 


(85938) 1999 DJ 4 


0.74 


0.8 


0.4 


0.17 


55 


15 


9.0 


2000 UGn 


0.77 


0.4 


0.2 


0.08 


280 


18 


22 


(3671) Dionysus 


1.2 


3.8 


1.5 


0.3 


190 


720 


120 


2002 CE 26 


0.67 


5.1 


3 


0.21 


88 


1300 


110 



Table 1: Ages of NEA binaries based on assuming their semimajor axes have evolved from 
much smaller initial values [eq. (1221) ]. Comparison for monolithic (k = 3/2/i), fluid (k = 3/2), 
an d rubble pile (k = 3/ 2 p, ru hhi. P ) strength for primary. Binary parameters from compilation 



bv lWalsh & Richardson! (12006af ). 



